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Abstract 



o 

(N 

Q I Oscillations of the F{R) dark energy around the phantom divide line, ujde = —1, both 

^ • during the matter era and also in the de Sitter epoch are investigated. The analysis during 

I the de Sitter epoch is revisited by expanding the modified equations of motion around 

\^ • the de Sitter solution. Then, during the matter epoch, the time dependence of the dark 

energy perturbations is discussed by using two different local expansions. For high values 
^ ■ of the red shift, the matter epoch is a stable point of the theory, giving the possibility to 

expand the i^(i?)-functions in terms of the dark energy perturbations. In the late-time 
^ i matter era, the realistic case is considered where dark energy tends to a constant. The 

{3JQ| results obtained are confirmed by precise numerical computation on a specific model of 

exponential gravity. A novel and very detailed discussion is provided on the critical points 
^ in the matter era and on the relation of the oscillations with possible singularities. 

00 ! 1 Introduction 



Modified gravity constitutes a very powerful and natural possibility to unify the physics of the 



0|0 I early time inflation epoch with that of the late-time acceleration stage, under the frame of a 

common theory (for a review, see |1]) and this without the need to introduce any extra fields 
(scalar, spinor, or other) as dark components. Among the possible versions of modified gravity 
there are in fact a number of them which comply without problems with all the local tests 
(for recent reviews, see OS]) and the most recent and accurate large-scale observational data 
(recent discussions can be found in [HEll^). There is, moreover, a serious expectation that 
5^ I all of the different stages in the universe evolution: inflation, radiation/matter dominance, 

the present and future acceleration epoch, will be qualitatively well described in frames of a 
specific class of modified gravity [1]. Furthermore, modified gravities can be reconstructed to 
give rise to other realistic scenarios of the universe evolution, such as the ekpyrotic one [7] 
and the so-called Little Rip cosmology [8j (for explicit examples of the late, see 

The present paper is devoted to the study of F{R) gravity properties in the matter and 
in the dark energy epochs. We will use the so-called fluid representation of F(R) gravity |10j . 
where the corresponding equations of motion are presented in the standard gravity, FRW form 
but with the addition of an effective fluid. Furthermore, this F(i?)-fluid will here contain an 
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explicit dependence on the gravitational terms. Specifically, we will study the oscillations that 
are generated in the equation of state (EoS) parameter wde around the phantom 

divide crossing, during the matter dominance era and in the de Sitter epoch (see also [llj). 
The evolution of the F{R) dark energy in these regimes will be here carefully investigated, 
both analytically and numerically. 

As is known, an important viable condition in modified gravity is the positivity of the 
second derivative of the cosmological function {F"(R) > 0) during the matter dominated 
era. This condition arises from the stability of matter perturbations in the weak-field regime, 
taking the curvature to be locally constant. When F"(R) < perturbations grow up and, as 
a consequence, the theory becomes strongly unstable |29| . In this paper a detailed analysis of 
such perturbations will be carried out, in the approach discussed, by considering the dynamical 
behavior of dark energy during the matter era. Since in the late-time matter era it is not stable, 
an expansion of the EOM with respect to the matter curvature is not allowed, and it will be 
necessary to introduce some physical assumptions in order to find correct results. Accurate 
viable conditions for the stability of the cosmological perturbations during the matter era 
will be found, in agreement with the constraint F"{R) > 0. We will see that an oscillatory 
behavior of the dark energy can occur during the matter era. This phenomenon leads to 
oscillations oi ujde, whose frequency and amplitude decay with the red-shift until the current 
accelerated epoch, while they substantially grow up as a consequence of perturbations around 
the de Sitter solution of the field equations. 

It will be proven, in particular, that the dark energy density approaches a stable de Sit- 
ter point, either as a positive power function of z + 1 or with an oscillatory behavior which 
amplitude decreases as {z + 1)^/^. In the late-time matter era, 1 < z < 3, a realistic, expo- 
nential gravity theory will be investigated numerically. We will actually obtain the precise 
time dependence of the dark energy perturbations. The analysis of this dark energy during 
the matter epoch will be performed using two different local expansions. For high values of 
the redshift the matter epoch will turn out to be a stable point of the theory and, thus, it will 
be possible to expand the F(i?)-functions in terms of the dark energy perturbations. In the 
late-time matter era, on the other hand, we will consider the realistic case where dark energy 
tends to a constant. The results obtained will be carefully checked by numerical computation 
on a specific model of exponential gravity. The corresponding attractors for this exponential 
model will be identified. 

We include a detailed discussion on the critical points in the matter era. In this set up 
only the perturbation associated with the scalaron will be investigated, anyhow the results 
obtained may be interesting and useful. They will be expressed both in semi-analytical and 
in numerical form, and we will prove their consistency with each other. We do not know of 
any other competing study on this quite non-trivial subject. As last point, some problems 
appearing in relation with the big values generated for the dark energy fluctuations, for large 
values of the red shift, will be discussed, by considering the relation between the oscillations 
obtained and possible future singularities. 

2 cjD£:-oscillations in realistic F{R) models of the DE- universe 

The action of modified F{R) theories with matter is [1] 




(1) 
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where g is the determinant of the metric tensor gi^j^y, 

£ (matter) jg ^j^g matter La grangian and 
F(R) a generic function of the Ricci scalar, R. Throughout the paper we wih use units where 
ks = c = h = 1 and denote the gravitational constant = SttGn = STr/Mp^, with the Planck 
mass being MpL = G^^'^ = 1.2 X lO^^GeV. We shah also write 

F{R) = R + f{R). (2) 

Thus, the modification to gravity is encoded in the function F{R), which is added to the 
classical term R of the Einstein-Hilbert action of General Relativity (GR). In what follows we 
will always discuss modified gravity in this form, by explicitly separating the contribution of 
GR from its modification. 

Taking the variation of Eq. ([1]) with respect to g^u, one has: 

- \r9,. = p^fT'''' + ^Tp^ [la^uiFiR] - RF'iR)] + (V^V, - 5^.n)F'(i?)| . 

(3) 

Here Rf^^^ is the Ricci tensor, the prime denoting the derivative with respect to the curvature 
scalar R, is the covariant derivative operator associated with the metric 5^,^, and = 
S'^'^V^Vjyc/) is the d'Alembertian of the scalar field (p. Finally, T^^^^^^ is the matter stress- 
energy tensor. We will work in the spatially-flat Friedman-Robertson- Walker (FRW) space- 
time described by the metric 

ds^ = -dt'^ + a\t)d^\ (4) 
where a{t) is the universe scale factor. The Ricci scalar is given by 

R = 6(2H^ + h) . (5) 



From ([3|) the following gravitational field equations are obtained: 

3F'iR)H^ = K^pm + \{F'{R)R - F{R)) - 3HF'{R) , (6) 

- 2F'{R)H = k\p^ + Pm) + F{R) - HF{R) . (7) 

Here, H = a{t)/a{t) is the Hubble parameter and the dot denotes time derivative d/dt, 
while p and p are, respectively, the matter energy-density and pressure. The following matter 
conservation law results 

p^ + 3H{pm+p) = Q. (8) 

For a perfect fluid, it gives 

p = u}p, (9) 

CO being the thermodynamical EoS-parameter for matter. For standard matter, u = and 
p = a{t)~^ while, for radiation, cj = 1/3 and p = a(i)~^. 

Eqs. ©-(IT]) can be written as functions of the effective energy density, Peff, and effective 
pressure, Peff, of the universe. This yields 

Peff = , (10) 

Pes = -^(2H + 3H^) , (11) 
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where 



Pcff 



F'{R) 



P + 



F'{R) 



P + 



2k2 



J-Uf'{R)R-F{R))-6HF'{R) 
-{F'{R)R - F{R)) + AHF'{R) + 2F'{R) 



(12) 



(13) 



We note that Eqs. (|10p -(|TT|) have in fact the same form as the Friedmann equations of General 
Relativity. Notice that the part corresponding to modified gravity has been formally included 
into peff and Peff- Hence, one must not forget that gravitational terms enter here in both 
sides of the FRW equation. 

We will be interested in the cosmological behavior of realistic models of modified gravity 
describing the de Sitter epoch of the universe today. The tag 'realistic' has been defined in 
previous papers and has to do with the feasibility of the models in view the all the most 
recent and accurate observational data. We define the dark energy density p£,E as p^E = 
Peff — p/F'(R) and introduce the variable [T2] 



Pde 



(0) 



^-(z + lf-xiz + l) 



(14) 



Here, p^ is the energy density of matter at present time, is the mass scale 



1^ Prn 



~ 1.5 X 10^^'^eV' 



and X is defined as 



where p^^ is the energy density of radiation at present, z the redshift parameter, z = l/a(t) — 1, 
and yH{z) is written as a function of z. 

The EoS-parameter for dark energy, lobe, is 



.(0) 



X 



Pm 



~ 3.1 X 10' 



(15) 



(16) 



WDE = -1 + -{z + 1 — ^ . 

3 yH[z) d{z) 

By combining Eq. (jlOp with Eq. ([5]) and using Eq. (|14p . one gets 

d^yniz) dyniz) i ^, r n 

+ Ji — T h J2yH{z) + Js = 0. 



dz 



where 



J- 



{z + l) 
1 



1 



1 - F\R) 



yH + {z + If + x{z + 1)^ &m^F"{R) 
1 2 - F\R) 



J3 = -3(z + l) 
Thus, we have 



[z + 1)2 \yH + {z + 1)3 + x{z + 1)4 2,m'^F"{R)^ 
(1 - F'(i?))((z + 1)3 + 2x{z + If) + {R- FiR))/i3m^) 



[z + lfiyH + iz + lf+xiz + m 



6rh'^F"{R) 



R = 3m^ [4yH{z)-{z + l) 



dyujz) 
dz 



+ {z + iy- 



(17) 

(18) 

(19) 
(20) 
(21) 

(22) 
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where we have used —iz + l)H{z)d/dz = H{t)d/d{lna{t)) = d/dt. 

In general, Eq. (jlSp can be solved in a numerical way, once we write the explicit form of 
the F(R)-model. Neglecting the contribution of matter, the trace of Eq. ([3]), 

3DF'{R) + RF'{R) - 2F{R) = , (23) 

gives us the well-known de Sitter condition 

2FiRds)- RdsF'{Rds) = 0, (24) 

where the Ricci scalar R^s is a constant. 

As a check, let us now study perturbations around the de Sitter solution of the dark energy 
density, to see that we are able to recover well known results. Here we restrict our analysis 
to homogeneous perturbations. The behavior of general, linear, inhomogeneous perturbations 
has been discussed in [13j, where the equivalence between the two approaches has been shown 
explicitly (see also the independent proof contained in [14j). 

The starting point will be 

yH{z)^yo+y{z), (25) 
where yo = R^g /12'm'^ is a constant and \y{z)\ « 1. Eq. (j22p leads to 

R = (^4y^ + 4y(^) _ + 1)^^^^ + + 1)') • (26) 

In this case, by neglecting the contribution of radiation and assuming the matter one to be 
much smaller than i/q, Eq. (fT8|) becomes 

'''^'^ + + ^A^yiz) = ^Siz + 1) , (27) 



dz^ (z + l) dz (z + l)2- 



where 



a = -2, (28) 
^" ^^RF"{R,s)' 

5-1 + 1—^^ (30) 

^~^RdsF"{Rds)- ^ ' 

Here we have performed the variation of Eq. (|26p with respect to R, and we have used 
Eq. (IMl) too. The solution of Eq. is 



,(.) = Co(. + i)K-"^^/(^=^) + |(. + 1)3 , (31) 

where Co is a constant. It is easy to see that \y{z)\ « 1 when z — — l"*", and, therefore, the 
de Sitter solution is stable provided that (see for instance [13 1 [16]) 

^'^^"''^ > 1 . (32) 



RdsF"{Rds) 
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We have two possible behaviors for the dark energy density in viable models of modified 
gravity, for a stable de Sitter universe. If 

25 F'{Rds) 

16 > R,sF"{R,s) ^ ^ ' ^^^^ 
the solution approaches the de Sitter point as a power function of (z+l), that is y(z) ~ (z+l)'^ , 
7 > 0. Otherwise, if 

F'{Rds) 25 
RdsF"{Rds) ^ 16 ' ^ ' 

the discriminant in the square root of Eq. (|3ip is negative and the dark energy density shows 
an oscillatory behavior whose amplitude decreases as (z + 1)^''^, when z — t- — 1^. As a 
consequence, we can write yH{z) as 

- 1^ ^ { p^iH.^) - H,sF"(R^) - ') " ^ ^ " ^ ^ 



A and i? being constants which depend on the boundary conditions. 
Using Eq. (fT7|) . we can evaluate the ojde parameter 

A5{z + lf 1 (z + l)T 
u:de = -1 + -t- - + 1^1- 36 

where 

1 



7 = - (1 - a ± V(l - «)' - 4/3j . (37) 



In the case of oscillating models which satisfy Eq. (|34p . one has 



with 



and 



= ( '^''»1f\-?lA^ (40) 



2 V \RjsF"(Rds) 1 

We observe that a;^)^; exhibits the same oscillation period of yniz) and that its amplitude is 
amplified by its frequency, written in the coefficients A' and B'. 

In the case of models of Eq. ([2]) , the oscillatory condition leads to 

1 + f'^^'^s) ^ 25 

Rdsf"{Rds) 16' ^ ' 

The large class of one-step models |17l I18| which mimic the cosmological constant in the high- 
curvature regime (that is, when f{Rds) — — 2A, where A is the cosmological constant) satisfy 
this condition {f"{Rds) << 1) and the corresponding uj^e oscillates around the value —1. 
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2.1 Time evolution 



Let us now consider the stable de Sitter solution of Eqs. ()25p and ()3ip in the case of 7 being 
a positive real number (what means that the condition (|33p is satisfied): 



yH{z)=yQ + Co{z + ir + j{z + lf . 



(42) 



Here (3, 5 and 7 have been given by Eqs. (I29D-(I3DD and yo = Hl^/rh^ = Rds/l^rh^ and 

Co is a constant. We will assume Cq > 0. The first equation of motion (|10p leads to 



= yniz) + {z + If = yo + Co(z + If + 



m 



1 



F'iRas) - RdsF"{Rds] 



iz + 1) 



(43) 



We will explicitly solve H as a function of the cosmic time t. By writing z + 1 as l/a(t), one 
gets 



a{t) 



1 



a{t) 



(44) 



We have here omitted the matter contribution to yniz). By considering t > 0, the general 
solution for the expanding universe is 



ait) 



Cpm^ 



sinh (^—^it + 



being (p a positive constant. It is then easy to obtain 



a{t) = a^e 



1 



(45) 



(46) 



where oq is a constant which depends on (j). As 7 > 0, we get a{t) ~ a^e^'^^^. For the Hubble 
parameter, one has 

^" (47) 



H = Hds coth [^-Hdslt + 

and, in general, for t > it is ~ H^s- 

Consider now the case of an oscillatory behavior followed by the condition (|34p . Then 

'a{ty 



H 



a{t) 



Hjg + m'^a{t)~2 (A cos (-Flog(a(t)-i)) + 5sin \og{a{t)-^))) , (48) 



where we have used Eq. ([35|) . omitting matter contributions, and the frequency T is 

F-- 



' ^F'jRds) 25 
RdsF"{Rds) 4 



Assuming a{t) ~ exp(i?o*)) Eq. ([l8]) yields 



H 



-A cos {FHQt) + B sin {FHot)] . 



Also, in this case, H ~ H^s- We should stress that, as 



■H{t) 



-H{z){z + r 



d 
dz 



H{z), 



(49) 



(50) 



(51) 



it is quite simple to see that any derivative of H{z) becomes singular around the zeros of the 
sinus and cosinus functions. 
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2.2 Singularities and final attractors 



We know, that if condition ([32]) is satisfied, the de Sitter point is stable, and it may be a final 
attractor of the system. However, in many papers, it has been demonstrated that the presence 
of finite-time singular solutions in expanding universe influences this behaviour. Finite-time 
singularities appear when Hubble parameter has the form 

where /iq is a positive constants and t < to. If /3 > — 1, or its derivative and therefore the 
Ricci scalar diverge and become singular at the time to- 

If a modified gravity theory shows such kind of singularity where R — t- ±00, the cosmologi- 
cal expansion could tend towards such asymptotic solution and the universe becomes singular. 
This is because singular solutions often are energetically accessible for the system and desta- 
bilize the model. As a qualitative example, we can consider the case of a realistic /(i?)-model, 
namely the Hu-Sawiki Model |17) . able to reproduce the de Sitter phase of our universe, 

^ ' C2{R/m^r + 1 " C2 C2{R/m^r + 1 ■ 

Here, is a mass scale, ci and C2 are positive parameters and n is a natural positive number. 
The model is very carefully constructed, such that c\rr? — 2A, where A is the Cosmological 
Constant, and in the high curvature region the physic of ACDM Model can be found. The 
Hu-Sawiki Model could become singular when R diverges and /3 = —n/{n + 2) > —1. 
Let us rewrite Eq. (j23p neglecting the contribute of matter in the well known form 



where 



dV^ 1 



aF'(R) -,l^nm-RF'iR)]. (55) 

Here, F'{R) is the so-called 'scalaron', this is the effective scalar degree of freedom of modified 
gravity, and VcS is the effective potential associated with the scalaron. In the case of Hu-Sawiki 
Model, the scalaron reads 



F'{R) = 1 ^'^^/^^ • (56) 



It tends to a constant when R — )■ ±00. Furthermore, by writing dV^s/dR as F"{R){dVcs/dF'{R)), 
in principle one can evaluate the potential VeS of Eq. (|54p through an integration. When 
R —7- ±00 one easily finds 



m^ci/c2 
?>C2{R/rn^ 



V,s{R ^ ±00) ^ ;„7-- 2,n (^ + 1) • (57) 



The Hu-Sawiki Model exhibits a stable de Sitter solution in vacuum, that may be the final 
attractor of the system. However, we can observe that if a singular solution with R diverging 
exists, it is at a finite value of V^s (in particular, it tends to zero) and the scalaron F'{R) can 
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crossover the potential in some point of cosmological evolution and arise the value F'(R) = 
for which catastrophic curvature singularity emerges. In general, it is possible to see that 
singularities violate the strong energy condition (SEC) describing acceleration. This is the 
reason for which realistic models of modified gravity describing the current acceleration of 
Universe could become unstable and fall into a singularity. The study of the singularities is 
fundamental in order to achieve a correct description of the Universe. As regard this argument, 
we will make a confront with realistic exponential model in the last section. 



3 a;£)£;- oscillations in the matter era 

Let us now consider matter era solutions. The trace of Eq. ([3]) reads, in this case, 

DF'{R) = ^ {2F{R) - RF'{R) - K^p^) , (58) 
where we have explicitly written the standard matter contribution to the stress energy tensor 

T-rmatter 

The critical points associated with matter dominated era for a generic F{R) model have 
been investigated in Ref. [20) . The results are 

F'{R) = 0, pr = 0, R = 3H'^ (59) 
F{R) = RF'{R). (60) 

As a consequence, from Eq. (|12p around the critical points, one has p^ff ~ Pmatter / F' [R) . 

Now we consider the perturbations of dark energy during the matter epoch, by assuming 
that yniz) <C (1 + z)^. Eq. ([TH]) reads, to first order in yH{z)/{z + 1)^, 

y^^'^ - (ITT) + (iw I RF"{R) ) = 

where R is written in full form, as in Eq. (j22p . and we have used the conditions (j60p . Here 
and in what follows we have put R — F{R) ~ 6m^y/f(z), as a consequence of Eq. ()12p and 
conditions (|60p . 

In the standard cosmological scenario, the efi^ects of dark energy are completely neglected 
when z ^ 3. In this case the Ricci scalar simply reads 

R = 3m'^iz + lf , (62) 

and Eq. ()6ip can be expanded as 

. . ^ , / . ^ 1 [7 {1-F'{R))F'"{R) \ 

. . 1 , 1 , 2(1 - F'{R))F"'{R) \ ( 2- F'jR) - F{R)/R \ 

^^^^^ {z + 1)2 V RF"{R) + F"(i?)2 )~V 2RF"{R) ^ + ^^ 
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In order to solve Eq. (|63p we can set z = zq + {z — zq), where \z — zo\ <^ zq, and perform a 
variation with respect to z. To first order in (z — zq), we find 

1 / 7 (1 - F'(fl„))F"'(fio! 



(zo + 1)2V RoF"{Ro) F"{Roy 
2-F'{Ro)-F{Ro)/R, \ 

^ + 2i?oF"(i?o) j ("° + ^) + 

/I 5F{Ro)/Ro - F'jRo) - 4 (2 - F(i?o) - F(i?o)/i?o)F^^(i?o) \ _ 
V2 QRoF"{Ro) 2F"{RqY J^^ 

where i^o = 3m^(zo + 1)^- The solution of this equation is 

yniz) = a + b-{z-zo) + C- e^("±^/^^)(^-^°) , (65) 

where C is constant and 

1 \ 6RlF"{Ro) + {2-F'{Ro))Ro-F{Ro 



Qrh^J 1 + 2RoF"{Ro) + 2(2 - F'{Ro) - F{Ro)/Ro)RoF"'{Ro)/F"{Ro) + ^ ^ 

/ Rl \ 7F"{Ro f + (2 - F^(i?o) - F{Ro)/Ro)F"'iRo) 

{4m^J [2RoF"{Ro)^ + F"(i?o) + 2i?o(2 - F'(i?o) - FiRo)/Ro)F"'iRoW 
[RF"{R^f + {bF{R^)/Ro - F'{R^) - 4)F"{R^)/3 - Ro{2 - F'{R^) - F{Ro)/Ro)F"'{Ro)] , 

Rq RoF"{Rof + {5F{Ro)/Ro - F'{Ro) - 4)F"{Ro)/3 - (2 - F'{Ro) - F{Ro))RoF"'{Ro)) 



2m^{zo + 1) 2RqF"{R^Y + F"{Ro) + 2(2 - F'{Ro) - F{Ro)/Ro)RoF'"iRo) 

(67) 

7 {l-F'{Ro))F"'{Ro) 

"=2 + 2FW ' ^ ^ 

1 2(1 - F'jRo)) F"'{Ro) 

^ - ^ + i?oi^"(iio) + ■ ^^^^ 

Let us now analyze this result. Since in the expanding universe (z — zq) < 0, it turns out that 
the matter solution is stable around Rq if a > and /3 > 0. This means that 

il-F'{R))F"'iRo) _7_ 

2F"{Ro)^ 2' ^ ' 

We can thus have an oscillatory behavior of the dark energy if the discriminant of the square 
root of Eq. ()65p is negative. 

3.1 Late-time matter era 

The effects of dark energy could be relevant at a late-time matter era, near the transition 
between the matter and de Sitter epochs (f ^ z <3). In this case, we cannot do an expansion 
of the i<'(i?)-functions in terms of yniz), as we did before. On the other hand, in realistic 
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models of modified gravity, yniz) tends to a constant value, as in Eq. (j25p . yniz) = yo + y{z), 
where yo — -R(is/127ri^ is related to the de Sitter solution and \y{z)\ <C yo (in this way, we 
can reproduce the correct dynamical evolution of the universe, as in the ACDM model). As a 
consequence, we can actually perform the variation of Eq. (j6ip with respect to y{z), to obtain 



y"{z) + y'{z) 



1 



iz + l) 



3 + 



(^ + 1) 
1 

2F"{R)R 



7 _ a - F'{R))F"'{R) 
2 ~ 



1 - F'{R)) + 



2F"{Rf 

{R-F{R)) \ 
R ) 



^ yo + y{z) 4F'{R) - 3 ^ ^^2) 



{z + 1)2 RF"{R) 



where 

R = 3m2 [(z + 1)=^ + Ayo] . (73) 

Also in this case, we can take z = zq + (z — zq), where \z — zq\ <^ zq, and doing the variation 
with respect to z, we find, up to first order in (z — zq), 



y"{z) + y'{z)- 



1 



{zo + 1) 



3 + 



ZQ + 1) 

1 



7 {1-F'{Ro))F"'{Rq) 
2 2F"{RqY 



+ 



1 



2F"{Ro)Ro 



1 1 - F'{Ro 



2 2F"[RQf 



-F"'{Ro) + 



1 - F'{Ro) + 
1 - F'{Ro) 



6F"{Ro)Ro 



Rq — F{Ro) 

Rq 

{z - zo) , 



yo + y{z) 



{zo + 1)2 RoF"iRo 



+ 



(74) 



where Rq = 3fh?{{zQ + 1)^ + 4yo)- In this expression we have used the condition ()60p and we 
have set Rq — F(Rq) ~ 6m?yo. Owing to the fact that yQ <^ Rq, the variations with respect 
to 2/0 cire at least of first order in yo/RQ. The solution of this equation is 



y 

where C is a constant, and 
Rq 



yo = a, 

(z) = b-{z-ZQ) + C- («±%^^) (--o) 



6m2(4F'(i?o) -3) 



(6F"(i?o)^o + 2 - F\Rq) - F{Rq)/Rq) + 



(75) 



(76) 



^{7F"{Rq)' + 2 - F'iRQ) - F{RQ)/RQ)F'"iRQ)) x 



(i?oF"(i?o)' 

, 3(zo + l)2(i?oi^"(i?r^2 



(1 - F'iRQ))RQF"'iRQ) + (1 - F'iRQ))F"iRQ)/3) , 

_ (1 _ F'iRQ))RQF"'iRQ) + (1 - F'(i?o))i^"(i?o)/3) 



2F"(i?o) 



_ 7 (l-F^(i?o))F^^^(i?o) 
" 2 2F"(i?o)2 



/3 



1 



RoF"{Rq) • 

The solution is stable around i?o if o; > and /3 > 0. This means that 

{1-F'{Rq))F"'{Rq) 



2F"{RqY 



7 

>-2' 



(77) 
(78) 

(79) 
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The oscillatory behavior of the dark energy occurs when the discriminant of the square root 
of Eq. ()75p is negative. 

We observe that the expression ()75p is more accurate than ([65|) . In general, if the conditions 
(I79])-(l80]) are satisfied for Rq = 3m^{zo + 1)^ + Um^yo, the conditions (I7I|)-(IZn]) wih be also 
satisfied, provided it is possible to use the approximation Rq = 3m^(zo + l)'^- In particular, 
since F'{R) ~ 1, the condition 

is sufficient to obtain stability during the matter era and an oscillating behavior of the dark 
energy. Such condition has to be satisfied by models which mimic the ACDM Model {F(R) ~ 
R-2A), for which F"{R) 0+. 

3.2 oj^j^;- parameter in the matter era 

We can use the solutions ()65p or (|75p and the expression (|17p to evaluate the parameter code 
around z = zq during the matter era. We get 



1 

6a 



b{zo + 1) + C"^^f "'^e^ («± 



(82) 



In the case of the oscillating models, for which — 4/3 < 0, this equation reads 



b , , e 2(^0+1) 

bJBE ^ -1 + —{zq + 1) H 

6a ba 



2(.o + l) + 2(zo + l) 



where A and B are constant, proportional to the period y^4/3 — a^/2(zo + 1). 



(83) 



4 Oscillating behavior of the dark energy in a exponential model 

In Refs. |17^ [T8l I21| several versions of viable modified gravity have been proposed, namely 
models, which are able to reproduce the current acceleration of the universe within the bounds 
set by the most recent observational data. They incorporate a vanishing (or fast decreasing) 
cosmological constant in the flat (i? — )• 0) limit, and exhibit a suitable, constant asymptotic 
behavior for large values of R. The simplest one was proposed in Ref. |18j and can be given 
under the form of Eq. ([2]), namely F{R) = R + f{R) , where 

/(i?) = -2A(1 -e-^/^°). (84) 

Here, A ~ 10~^^eV^ is the cosmological constant and i?o ~ A a curvature parameter. In flat 
space, one has /(O) = and recovers the Minkowski solution. For R ^ Rq, f{R) — — 2A, and 
the theory mimics the ACDM model. The viability of this kind of A models has been studied 
in Ref. [13 [22]. We are allowed to set 

RQ = bA, (85) 
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where 6 > needs to be sufficiently small, in order to find the de Sitter solution Ri^s at high 
curvature {f{Rds) ~ ~2A), so that de Sitter condition (|24p leads 



Rds ~ 4A , 

in accordance with the ACDM model. The stability condition (I32p yields 



6 > 4. 



(86) 



(87) 



It is easy too see that this condition is always satisfied and, in particular, the oscillating 
condition of Eq. (jlT]) is satisfied for any value of b. In what follows, we will consider the case 
6 = 1. 

Eq. (jlSp can be solved in a numerical way. To this purpose we have imposed the following 
initial conditions at z = z^: 

dyn 



d{z) 



VH 



0. 



A 



which correspond to the ones of the ACDM model. This choice obeys to the fact that in the 
high redshift regime the exponential model is quite close to the ACDM Model. The value of Zi 
has been chosen so that RF"{z = Zi) ~ 10~^, assuming R = 3rh'^{z + 1)^ + 4A, and Zi = 2.6. 
In setting the parameters we have used the last published results of the VFMAP, BAO, and 
SN surveys [23j and also that A = {7.93)m^. 

A plot of yniz) is shown in Fig. 1. In order to study the behavior of dark energy, we also 
extrapolate the value of the corresponding density parameter Qde{z), as 



Pde 



yniz) 



Pes yH{z) + {z + lf + x{z + l)^' 
In Fig. 2, the corresponding plot of ilDE(-2) is overlapped with the one for yni^) — A/3m^. 



(89) 




Figure 1: Plot of yniz)- 
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Figure 2: Plot of Qde(-2) overlapped with yniz) — A/3m?. 

We must point out that the data are in perfect accordance with the most recent, accurate 
observations of our current universe, where 

n^EiO) = 0.721 ±0.015. 

In fact, this value is to be compared with what we get at present red shift, z = 0, namely 
Qde(O) = 0.732. 



4.1 De Sitter oscillations 

Dark energy is dominant near the de Sitter solution describing the present acceleration of the 
universe (for z < 1.5, Qde > 0.5). As the de Sitter solution is a stable point of the theory, it 
turns out that yniz — ?• —1) ~ Rds/^'^rh?, where R^s is the solution of Eq. ([24|) . In our case 
this equation is transcendental and such that taking f{Rds) — /(4A), one has 

Rds ^'iAx (0.947) . (90) 

A computational evaluation (extrapolated from yniz = 0.9)) gives R^s ~ 3.725A so that we 
are remarkably close to the corresponding value for the ACDM model, namely R^s = 4A. 
Further, as the condition ()34p is satisfied, we can predict an oscillatory behavior of yniz), 
which is shown in fact in Fig. 1. The value oi loqe oscillates infinitely often around the line 
of phantom divide at ojde = — Ij as in Eq. (l38|) . In Fig. 3 the behavior oi ijj£,E is depicted. 
In order to appreciate the constant frequency of code with respect to the logarithmic scale 
during the de Sitter phase, in Fig. 4 we plot the values of 



u;bE{log{z + 1)) 



12m2 / 1 



{Z + 1)~2 



(91) 

as a function of log(z-|- 1), for —0.999 < z < (here, R^s = 3.725A). In this way, we stress the 
oscillating part of Eq. (|38|) . whose frequency is proportional to 2it/ y^4/(Rdsf"iRds)) ~ 25/4 ~ 
1.570. As last point, we should remark that the amplitude oibJ^E is amplified with respect to 
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yniz) by its frequency, and it decreases as {z + 1)^/^. These data are also in good accordance 
with the observations of our present universe, which yield 

UDE = -0.97210:061 . (92) 
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Figure 3: Plot of 


ojde for — 1 < z < 2. 



(jj 
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Figure 4: Plot of ojde for log[-0.999 + 1] < log(z + 1) < log[0 + 1]. 



4.2 Matter oscillations 

In the high curvature region matter is dominant {z > 1.5 with Qde{z) ~ 0.1) and yu ~ 
A/3m^, as it is clear in Fig. 2. Eq. ()75p yields an estimation of the dark energy density value 
around zq. Owing to the fact that, in this region, f"{R) is very close to 0"*", dark energy 
oscillates as 

yniz) = a + e 2(^0+1) 
where A and B are constant, and 

a ~ A/Sm^ , (94) 



A sin 



[zo + 1) 



{z — zq) \ + B cos 



VP 



[zo + 1 



-{z - Zo) 



(93) 
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a 



3, 



A 



3m^((zO + l)3+4A/3m^) 



(95) 
(96) 



6m2[(^Q + 1)3 + 4A/3m2] 
Here we have used Eqs. d76|) -(f78]). We must stress that, by using our formula, it turns out 



that h ~ exp(— i?/A) ^ {yo/Rof', which we have not been able to evaluate, its value being 
extremely close to zero. The frequency J-{zq) of dark energy oscillations is 



A 



6m2[(zo + 1)3 +4A/3m2" 



3m2((zQ + l)3+4A/3A2) 1 

e A 



(^0 + 1) 



(97) 



while the amplitude decreases as exp(3/2(zo + 1))(^ — zq). We can verify the validity of our 
formula by analyzing in detail the graphics of yDE{z) in the vicinity of zq = 2.60, 2.55 and 
2.50. Such graphics are shown in Figs. 5-7. We have chosen an interval of \z — zq\ = 0.02. 
The period T{zq) of dark energy oscillations has to be T{zq) = 2'k/T[zq) and the number of 
crests in our interval is n\z~zo\{zQ) = 0.02/T(2;o)- The predicted values are 

r(2.60) ~ 0.003 , no.o2(2.60) ~ 6.667 ; 



r(2.55) ~ 0.004, 
r(2.50) ~ 0.006, 



no.o2(2.55) - 5.000; 
rio.02(2.50) ~ 3.333. 



These values are in good accordance with the numerical computation. The dark energy density 
IS very close to A/3m2 ~ 2.64333 and one can check that the amplitude of oscillation decreases 
with z. 




2.64333 



Figure 5: Plot of yniz) in the vicinity of zq = 2.60. 

Since the dependence of the amplitude J-{zq) on zq is weaker when zq < 2.5, we can study 
the oscillations of the dark energy density around zq = 2.45 and zq = 2.40 by using a longer 
interval |z — Zo| = 0.05. We find 

r(2.45) ~ 0.008 , no.o5(2.45) ~ 6.250 ; 

r(2.40) ~ 0.010 , no.o5(2.40) ~ 5.000 . 
These results can be compared with Figs. 8-9. 
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Figure 6: Plot of yniz) in the vicinity of zq = 2.55. 
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Figure 7: Plot of yniz) in the vicinity of zq = 2.50. 
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Figure 8: Plot of yniz) in the vicinity of zq = 2.45. 



When z < 2.40 the periods become too large with respect to the change of J-{zq) and we 
can no more distinguish oscillations in the dark energy until the beginning of the de Sitter 
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Figure 9: Plot of yniz) in the vicinity of zq = 2.40. 



epoch, when the dark energy behavior is governed by Eq. ()35p . On the other hand, the 
effects of such oscillations are amplified in the expression (j83p for uj^e, where the amplitude 
of oscillations is proportional to the period. Fig. 10 shows the behavior oi ujde inside the 
region 2.3 < z < 2.6. Since the frequency oi ujoe is the same as for yE)E{z), we can observe 
there how it decreases with z. 




Figure 10: Plot of ujde as a function of 2.3 < z < 2.6. 



5 Oscillations and singularities 

A finite-time, future singularity occurs when the Hubble parameter or some of its derivatives 
(and, therefore, the energy density, scale factor or, more simply, the Riemann tensor compo- 
nents) diverge |24j . The simplest well-known example of this future singularity is phantom- 
induced Big Rip [27]. The softer kind of such singularities are of Type IV. A Type IV singu- 
larity appears if the Hubble parameter and the Ricci scalar R (which depends on H and H) do 
not diverge, but higher derivatives of H do tend to infinity. Generally speaking, singularities 
are called 'finite-time' if they occur in a finite time during the cosmological evolution of the 
universe, yielding nonphysical components of the Riemann Tensor. 
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Let us return to exponential gravity models. By fitting parameters, they could become 
indistinguishable with respect to the ACDM model. They mimic, in fact, the cosmological 
constant to high precision and ujde can be made extremely close to the value of —1. The 
transition crossing the phantom divide does not cause any serious problem to the accuracy of 
the cosmological evolution arising from such kind of models. In Fig. 11 we depict the behavior 
of the Ricci scalar for the exponential model of the previous chapter. We can see in fact how 
it decreases with the red-shift, as one would expect, and the amplitude oscillations of the dark 
energy density are so small that one cannot observe them. 



R/A 




-0.5 0.5 1.0 1.5 2.0 2.5 



Figure 11: Plot of i?/A as a function of z. 

On the other hand, we have seen how, in a model of this kind, f"{R — )• oo) — )• 0"''. As a 
consequence, the frequency J-{zq) of the dark energy oscillations increases with the redshift, 
as in Eq. ([97|) . Since the method we have used to analytically evaluate such oscillations yields 
a linear expansion of the EoM, we cannot observe any singularity in the time derivative of 
yniz) and, correspondingly, in the Hubble parameter. But we have that 

oc iHzoT , (98) 

where to is the time at redshift zq. For extremely large values of the redshift, the high 
time derivatives of the Hubble parameter become infinitely large and approach an effective 
future singularity, showing a different behavior of the exponential model with respect to the 
ACDM one. This fact suggests that, back into the past, other terms could have become 
important in modified gravity theories. For example, in Ref. |25| [26] it is shown how terms 
proportional to protect the theory against future time singularities and yield a finite value 
for f"(R). Simultaneously, these terms induce also an early-time inflation, thus unifying 
primordial inflation with dark energy i28j 

A simple modification of the exponential model which incorporates the inflationary era is 
given by a combination of the function discussed above with another function, which repro- 
duces the cosmological constant during inflation. A natural possibility is [19J 

F{R) = R-2a(^1- e~^) - A^ M - e"("ft) j + 7i?" , (99) 
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where Ri and Aj assume the typical values of the curvature and expected cosmological constant 
during inflation, namely Ri, Aj ~ 10^'^~^^eV^, while n is a natural number larger than one. 
The last term of such equation jR", where 7 is a positive dimensional constant and a a real 
number, is actually necessary in order to obtain an exit from inflation. Typically, 7 ~ l/i?°~^ 
and a > 2, so that the de Sitter solution of inflation results to be unstable and the effects of this 
term vanish during the matter era, when R<^ R^. The detailed analysis of this extended class 
of exponential models for early- and late-time acceleration has been carried out in Ref. [19] . 
Here we restrict ourselves to note the effects of the power-like function on the dark energy 
fluctuations during the matter epoch. Now, Eq. (|97p reads 



' ' (100) 



where Rq is the curvature at redshift zq. Since a > 2, while the exponential term becomes 
small on approaching the effective singularity aX, Rq ^ 00, the power- like term becomes larger, 
and the frequency of dark energy oscillations does not diverge. 



6 Attractors of exponential models 

The exponential model (j99p has three solutions (critical points) of the theory in vacuum: 
an unstable de Sitter solution for R^s ~ Ri (for more detailed evaluation see Ref. |19| ) 
describing inflation, the second one, which is stable, at Rds — 4A, and the last one, at i? = 0, 
corresponds to the Minkowski solution of Special Relativity. When, in early time, the curvature 
of the universe is of the same order of the system gives rise to a de Sitter solution where 
the universe expands in an accelerating way but, suddenly, it exits from inflation and tends 
towards the minimal attractor at Rds = 4A, unless the theory develops a singular solution for 
ii —)• 00. In such case, the model could exit from inflation and move into the wrong direction, 
where the curvature would grow up and diverge, and a singularity would appear. It has been 
demonstrated that exponential gravity is actually free from singularities. However, in the very 
asymptotic limit, when R^ Ri, the model reads 

F{R > Ri) ~ 7i?" . (101) 

A class of asymptotic solutions of Eqs. ([Gj-dT]) in vacuum, for the model of Eq. (jlOip . at the 
limit t — )• 0"*" results 

H{t) = ^ , (102) 

with ^ 

R^12^. (103) 

Here, Hq is a large positive constant and /? a positive parameter so that /3 = 1 or /3 > 1. 

This result shows that in the limit R — )• -|-oo the model exhibits a past singularity, which 
could be identified with the Big Bang one. However, with a change t — )• {Iq — t), we may 
obtain a future-time singularity solution. It is important to stress that this kind of solution 
is disconnected from the de Sitter phase of inflation, where the term R is of the same order of 
■jR", and is therefore not negligible as in Eq. (jlOip . In the very asymptotic limit, the scalaron 
F'{R) results 

F'{R) = -faR'^-^ . (104) 
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We can also evaluate the potential V^s of Eq. (fM|) through integration of F" {R){dVcQ / dF' {R)) = 
dVcfi/dR. One easily finds 

V;„(fl>>ft) = ^'°'3°-;l<^-°' fi-°-'. (105) 

We observe that, in order to reach the singularity, the scalaron has to cross over an infinite 
potential barrier {Ves{R — )• oo) — )• oo) and go to infinity iF'{R — )• oo) — )• oo). But clearly this 
dynamical behavior is forbidden. 

We may safely assume that, just after the Big Bang, a Planck epoch takes over where 
physics is not described by GR and where quantum gravity effects are dominant. When the 
universe exits from the Planck epoch, its curvature is bound to be the characteristic curvature 
of inflation and the de Sitter solution takes over. When inflation ends, the model moves to 
the attractor de Sitter point. In this way, the small curvature regime arises, the first term of 
Eq. ()99p becomes dominant and the physics of the ACDM model — as we have seen in the 
preceding sections — is reproduced. 

7 Discussion 

It is well-know that an important viable condition in modified gravity is the positivity of 
the second derivative of the cosmological function {F"{R) > 0) during the matter dominated 
era. This condition arises from the stability of matter perturbations in the weak-field regime, 
taking the curvature to be locally constant. If F"[R) < 0, such perturbations grow up and, as 
a consequence, the theory becomes strongly unstable [29] . In our approach, a detailed analysis 
of these perturbations has been carried out, by considering the dynamical behavior of dark 
energy during the matter era. Accurate viable conditions for the stability of the cosmological 
perturbations during the matter era have been found, namely Eqs. (|79p -(f80|). which are in 
perfect agreement with the constraint F"{R) > 0. We have also seen how an oscillatory 
behavior of the dark energy can occur during the matter era. This phenomenon leads to 
oscillations oi uoE, whose frequency and amplitude decay with the red-shift until the current 
accelerated epoch, while they substantially grow up as a consequence of perturbations around 
the de Sitter solution of the field equations. 

More specifically, we have investigated the oscillations of dark energy around the phantom 
divide line, ojde = — Ij both during the matter era and also well in the de Sitter epoch. The 
analysis of the dark energy during the de Sitter epoch could be carried out by expanding the 
modified equations of motion around the de Sitter solution. It has been seen that the dark 
energy density approaches a stable de Sitter point either as a positive power function of z + 1 
or through an oscillatory behavior whose amplitude decreases as (z + 1)'^/^. We have also 
found the time dependence of the dark energy perturbations. The analysis of the dark energy 
during the matter epoch has been performed by using two different local expansions. For high 
values of the red shift, z ^ 3, the matter epoch is a stable point of the theory and, thus, it is 
indeed possible to expand the i^(i?)-functions in terms of the dark energy perturbations. On 
the other hand, in the late-time matter era, 1 ^ ^ 3, we have considered the realistic case 
where the dark energy tends to a constant. Then, the results obtained have been accurately 
confirmed by a numerical computation on a specific model of exponential gravity. 

As last point, some problems appearing in relation with the big values generated for the 
dark energy fluctuations, for large values of the red shift, have been discussed in the paper in 
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detail. When F"{R) is very close to zero, simultaneously with the fact that modified gravity 
mimics very well the ACDM model, these fluctuations approach an effective singularity in 
high derivatives of the Hubble parameter, which become infinitely large. This fact suggests 
that a correction of the Einstein equations in the small curvature region is related with a 
modification of gravity at high curvature, produced by the leading terms of inflation. Our 
result can be easily generalized to other different models of dark energy, as the Little Rip one, 
where the effective wde approaches the value —1 for i — )• oo asymptotically. Note that our 
approach can be applied to other modified gravities, too (see, for instance, the review [30j). 

We have provided a detailed discussion on the critical points in the matter era. In this 
set up only the perturbation associated with the scalaron have been investigated, anyhow the 
results obtained are quite interesting and useful. They have been expressed both in semi- 
analytical and in numerical form, and we have shown in the paper that they are strictly 
consistent with each other. 
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